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A short personal account is given of the impact of HERA data and the 
influence of Jan Kwiecinski on low x physics. 

1. Introduction 

That HERA data was to change our understanding of pQCD was an- 
nounced at the Durham Phenomenology Workshop on HERA Physics in 
1993 (the first of three devoted to HERA Physics, all with a large input 
from JK). Albert De Roeck presented the first preliminary HI data from 
HERA on F 2 measurements at x < 0.01, below the fixed target region. 
Within large errors that data showed a rising F 2 as x decreased - the qual- 
ity of the data may be judged from Fig ^ showing comparable ZEUS data. 
The result was the major talking point of the 1993 HERA workshop and 
discussion has continued almost unabated since then. One gets a sense of 
why the result was surprising by comparing Fignwith FigEJ The latter plot 
shows a range of predictions from a paper by Askew, Martin, Kwiecinski & 
Sutton (AKMS) [2 with F 2 data from the NMC and BCDMS experiments 
for x > 0.01. The least model-dependent extrapolation from the measured 
data would appear give a 'flat' F2 as x — > as shown by the dash-dotted 
curve, but what do the rising curves represent? 

2. F 2 in pQCD 

A rising F 2 at small- x is predicted by pQCD, but the predictions do not 
give a scale in either x or Q 2 at which one might expect to see such effects. 
In the DGLAP 1 formalism at leading order the gluon splitting functions 

P P (1) 



1 Dokshitzer, Gribov, Lipatov, Altarelli & Parisi, collinear factorization. 
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Fig. 1. The first ZEUS F 2 data from HERA, published 1993 PJ. 

are singular as z — > 0. Thus the gluon distribution will become large as 
x — > 0, and its contribution to the evolution of the parton distribution 
becomes dominant. In particular the gluon will 'drive' the quark singlet 
distribution, and hence the structure function F%, to become large as well, 
the rise increasing in steepness as Q 2 increases. Quantitatively, 



dg(x,Q 2 ) ^ a s (Q 2 
dlnQ 2 



2tt 



/ — g{y,Q 
Jx y z 



(2) 



may be solved subject to the nature of the boundary function xg{x,Q 2 ) ). 
Inputting a non-singular gluon at Qq , the solution is |IJ [3] 



xg(x, Q 2 ) ~ exp 2 



X 



(3) 



where 



(4) 
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Fig. 2. Calculations from AKMS J5] of the behaviour of F2 and Fl at low x based 
on the BFKL equation. The difference between the upper and lower continuous 
curves is in the cut-off imposed to control diffusion in transverse momentum {kq = 
1 , 2 GeV 2 for the upper, lower curves respectively. The dashed curves show the 
effect of including shadowing effects with a proton radius of R = 5, 2 GeV -1 . The 
almost flat dash-dotted curves are the contributions excluding the BFKL effects. 

Given a long enough evolution length from Qq to Q 2 , this will generate a 
steeply rising gluon distribution at small x, starting from a flattish behaviour 
of xg(x, Q 2 ) at Q 2 = Qq. 

Over the x, Q 2 range of HERA data this solution mimics a power be- 
haviour, xg(x,Q 2 ) ~ x~ Afl , with 

_/12 1n(t/t ) y 

where t = ln(Q 2 /A 2 ), to = ln(Qg/A 2 ). This steep behaviour of the gluon 
generates a similarly steep behaviour of F2 at small x, F2 ~ where 
A = X g - e. 
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Fig. 3. The BFKL gluon ladder diagram for deep inelastic scattering with a quark 
box at the top connecting to the 7* and the proton at the bottom. 

So why was the observed rise of Fi unexpected? Because, before the ad- 
vent of HERA data, the starting scale for perturbative evolution had always 
been taken to be Qo~4GeV 2 , to be sure that a perturbative calculation 
would be valid. The evolution length from Q 2 of 4 to 15GeV 2 is not large 
enough to generate a steep slope from a flat input. 

However Jan Kwiecinski together with Martin, Roberts and Stirling |3] 
had already been considering alternatives to 'conventional' DGLAP with flat 
input distribution. The KMRSBo and KMRSB_ parton distribution func- 
tions of 1990 were both compatible with the data: KMRSBo had a con- 
ventional flat gluon distribution input at Qq = 4GeV 2 ; whereas KMRSB_ 
had a singular gluon distribution x -0 ' 5 , even at low Q 2 (= Ql). 

Why should one consider a steep input gluon?. Because at small x terms 
in ln(l/x) are becoming large and the conventional leading In Q 2 summation 
of the DGLAP equations does not account for this. It may also be necessary 
to sum leading ln(l/x) terms. Such a summation is performed by the BFKL 2 
equation. To leading order in ln(l/x) with fixed a s , this predicts a steep 
power law behaviour 

xg(x,Q 2 )~f(Q 2 )x- x e (6) 



2 Balitsky, Fadin, Kuraev & Lipatov, fcr factorization. 
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where 

A„ = ^41n2~0.5 (7) 

(for a s ~ 0.2, as appropriate for Q 2 ~ 4GeV 2 ). 

These ideas were taken further in the paper by AKMS, see Fig. |^1 in 
which BFKL evolution was incorporated directly into the calculation of the 
parton densities. The rising curves show different calculations of F2 and 
Fl from the BFKL equation, the dashed curves show the effect of including 
some damping by shadowing which will be discussed later. The dash-dotted 
curve is the prediction for the structure functions without BFKL effects. 

Thus the observation of a steep behaviour in F2 at a relatively low Q 2 was 
seized upon by those who saw this as evidence for BFKL. On the other hand 
it was quickly realised that by lowering the starting scale in the DGLAP 
approach to Qq rs 1 GeV 2 one could also fit the data without a singular 
gluon input distribution. This later approach owed much to the work of 
Gliick, Reya and Vogt (GRV) [B|. Higher statistics data taken during 1993 
showed that an effective power X g « 0.5 was actually too steep. 

Kwiecinski was involved with many of the approaches to find a credible 
phenomenology of the BFKL equation. One of the most significant of these 
is the need for a 'kinematic constraint' to control the diffusion of transverse 
momentum down the 'gluon ladder' in the BFKL approach (see Fig. |3J). 
Unlike the DGLAP summation in which the /ct of the ladder are strictly 
ordered, those in the BFKL summation are not and can diffuse downwards 
into the non-perturbative region. Consider a link in the gluon chain where 
the longitudinal momentum fraction decreases from x/z to x and the trans- 
verse momentum k' T changes to kx, with emission of a gluon of transverse 
momentum qx- One requires 

k 2 T jz > q\ (8) 

in order that the virtuality of exchanged gluons is controlled by their trans- 
verse momenta. This implies that k\jz > k'^, for any given value of kx- 
If one considers the effect of this constraint on the solutions for the BFKL 
equation one finds that it modifies the asymptotic solution x~ Xg , such that 
X g is reduced from X g ~ 0.5 to X g ~ 0.3. These ideas were built into an am- 
bitious approach by Kwiecinski, Martin and Stasto [Sj to combine the BFKL 
and DGLAP approaches in a single model to be applied to the HERA and 
fixed target data. The key idea was to use the BFKL kernel, with the kine- 
matic constraint, applied to the unintegrated gluon density (i.e. differential 
in kx as well as x) to give an input gluon density to DGLAP evolution. This 
gives rise to a pair of coupled equations for the gluon and qq sea which can 
be solved numerically. The small valence quark contribution was taken from 
the leading order GRV set 9 to give a complete model for F2, F^ harm and 
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Fig. 4. The combined BFKL-DGLAP model of KMS fit to F 2 data from ZEUS, 
E665, BCDMS & NMC. From 0. 



Fl- Only two parameters were needed to describe xg(x, fcg) at the starting 
scale. A very good fit to the HERA 1994 F2 data together with data from 
E665, BCDMS and NMC was obtained, as shown in Fig.H 

During the period 1994 - 1996 the conventional DGLAP NLO QCD 
global fits were also refined and gave very good x 2 fits to the HERA data, 
starting with a conventional flat gluon input, provided that a low starting 
scale Ql = 1 GeV 2 was used. In 1996 the HERA data from the 1994 run 
were published, showing that that the steep slope of F2 extended as far 
down in Q 2 as Q 2 ~ 1.5 GeV 2 . The precision of the data had now increased 
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sufficiently that the low-x behaviour of the sea and the gluon distributions 
could be fitted separately. These results led to a new kind of surprise since 
it became clear that the flat input which had been used for both the sea 
and the gluon (Ag = \ g ~ 0) was actually a compromise between a sea 
distribution which remains steep, As ~ 0.2 down to Q 2 ~ 1 GeV 2 , and a 
gluon distribution which is becoming valence- like, X g < at low Q 2 . Fig. [5] 
shows a recent result from a ZEUS global fit using the complete HERA-I 
data sample [7j. This behaviour contradicts the original argument that the 
steep behaviour of the gluon distribution is driving the steep behaviour of 
the sea which is measured in F2. 

These conclusion were strengthened when the data from the 1996/7 runs 
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were published in 2000/2001. It became possible to make detailed analyses 
of the experimental uncertainties on the parton distributions and thus put 
limits on the kinematic range of applicability of the DGLAP formalism. For 
Q 2 < lGeV 2 the gluon distribution becomes negative at small x, see Fig. \5\ 
and the fit can no longer describe the data, see Fig. H3 Furthermore, even 
though the DGLAP formalism works for Q 2 > 1 GeV 2 , doubts remain about 
its applicability when X g < A<j. With only one observable measured with 
high precision and all models depending on parameters fit to data, one has 
sufficient flexibility to get good descriptions in both the pure DGLAP and 
'BFKL enhanced' approaches. The striking rise in F<i does begin to abate 
for Q 2 below 2 GeV 2 as Fig. El illustrates. At this point is useful to bring in 
data on the slopes or derivatives of F2. 
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Fig. 7. The slope - * as a function of a; from ZEUS 1994/5 HERA data. From 
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3. F2 slopes 

The behaviour of F2 may be characterized in terms of two slopes, with 
respect to InQ 2 and m(l/x) respectively. The latter has already been used 

OF 

in the discussion of the gluon density at low x, here A(Q 2 ) = — — ; — -— or 

om(l/x) 

F2(x,Q 2 ) ~ x~ x (Q '. At low x the form of the DGLAP equations is such 
that at LO one has (very roughly) 

F 2 (x,QV*£(x,Q 2 ), and ~ P qg ® xg(x,Q 2 ). 

Thus the behaviour of di^/c? In Q 2 may come from either the gluon 
density or the splitting function or both. If it could be shown that any un- 
conventional behaviour was attributable to P qg then this would indicate the 
need for an alternative evolution equation. In any case it is well worthwhile 
the experimental groups calculating the two slopes with a proper treatment 

of systematic errors. — — —7 - the 'Caldwell plot' 3 as a function of x is 
omQ z 

3 A preliminary version of this plot was first shown by Allen Caldwell at the 1997 
DESY Theory Workshop. 
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shown in Fig. This is another result from HERA that caused a lot of 
interest and discussion! The plot shows the slopes calculated from the Fi 
shown in Fig. |1J An important detail to note is that because of the strong 
correlation between x and Q 2 inherent in deep inelastic kinematics, the 
slopes are not evaluated at a fixed values of Q 2 , but rather at the mean 
values shown along the top of the plot for each data point. The plot shows 
that the rate of rise of F2 as a function of Q 2 is starting to abate as x < 10~ 3 
and Q 2 < 8 GeV 2 . To understand what the plot might indicate and why 
one is also very interested in where and how Fi stops rising at low x it is 
helpful to look at the HERA data from another point of view. 

4. cr(7*p) and saturation 

The strong rise of F2 at low x was unexpected for another reason. At 
small x, F2 is related to the cross-section for 7*p scattering by 

a^(W 2 ,Q 2 )^^F 2 (x,Q 2 ), with W 2 « Q 2 /x, (9) 

where W is the ~f*p centre of mass energy. This relation implies that a 
rise of F2 at small- x corresponds to a rising 7*p cross-section with W 2 . 
The relationship above suggests another approach to understanding the be- 
haviour of F2 at low x and that is the framework of Regge theory used to 
explain the high energy behaviour of hadronic total cross-sections. Regge 
theory predicts that the high energy behaviour of hadronic scattering am- 
plitudes is lmA(ab — » cd) ~ YJ/3jS ai , where a, is the intercept of a Regge 

i 

trajectory which has the right quantum numbers for an exchange in the 
crossed channel ac — > bd. Using the optical theorem, a total cross-section 
will vary as a tot (ab) ~ For the corresponding forward elastic 

i 

scattering amplitude, the leading Regge trajectory has the quantum num- 
bers of the vacuum and is known as the Pomeron, for which the value 
ap ~ 1.08 has been determined from hadron-hadron data. This prediction 
also describes the high energy behaviour of photoproduction cross-section 
measurements successfully. The model was extended to describe virtual- 
photon proton scattering by Donnachie and Landshoff who assumed 
that the Q 2 dependence would reside in the residue functions Pi(Q 2 ) only 
and that the intercepts, oq would be independent of Q 2 . Through Eq. |§] 
this approach then predicts a flattish input for the gluon and sea PDFs, 
since (j(-y*p) ~ implies F2 ~ x x ~ ai = x -0 ' 08 . While the DL ap- 

proach successfully described the pre-HERA low x data for Q 2 values up 
to about 10 GeV 2 , it cannot describe the steeply rising F2 data measured 
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at HERA. The problem is summarised graphically in Fig. |H] which shows 
d In F 

HQ 2 ) = o-, , i x or ^2 ~ x~ x( - Q \ at low x. For Q 2 values less than 
ciln(l/i) 

1 GeV 2 or so, the value of A is consistent with that from hadronic Regge 
theory, whereas for Q 2 > 1 GeV 2 the slope rises steadily to reach a value 
greater than 0.3 by Q 2 ~ 100 GeV 2 . This larger value of A is not so far from 
that expected using BFKL ideas. Indeed the BFKL equation can be viewed 
as a method for calculating the hard or perturbative Pomeron trajectory 
- in contrast to the soft or non-perturbative Pomeron of hadronic physics 
with intercept around 1.08. Donnachie and Landshoff |12j have extended 
their Regge model by the addition of a hard Pomeron with intercept 1.44, 
which allows them to describe the low x HERA data up to Q 2 values of a 
few hundred GeV 2 . Theirs is one of many attempts, using ideas from Regge 
theory, to model the 'transition region' from real photoproduction through 
to deep inelastic scattering with Q 2 >> lGeV 2 . 

Either considering the rise of F<i as x — > or the steep energy dependence 
of <r(7*p) at large Q 2 , such behaviour will eventually violate the Froissart 
bound. If the origin of the rise arises from a high gluon density then this 
problem may be avoided because the gluons can shadow each other from 
the 7*. At even higher densities the gluons will 'recombine' via the process 
gg — > g - the inverse of gluon splitting - and the gluon density and hence F2 
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Fig. 9. Approaches to physics at low- a;. Courtesy of A D Martin 



will saturate. These ideas have been formalized by Gribov, Levin & Ryskin 
by the addition of a non-linear term to the equation for gluon evolution 



d 2 xg(x, Q 2 ) 
dlnQ 2 dln(l/x) 



—xg(x,Q 2 ) 

7T 



81a* 



WQ 2 R 2 



xg(x,Q 2 ) 



(10) 



When xg(x,Q 2 ) ~ irQ 2 R? / a s (Q 2 ) the non- linear term cancels the linear 
term and evolution stops, this is saturation. Shadowing effects were included 
in the calculations by AKMS shown in Fig. |2] and in the KMRS parton 
densities. 

The various evolution equations applicable across the x, Q 2 plane are 
summarised schematically in Fig. H 4 Note the appearance of the 'critical 
line' above which the gluon density is high enough for non-linear effects to 
be important. In this region and for large enough Q 2 one has the possi- 
bility that a s will be small enough for weak-coupling but non-perturbative 
methods to be applicable, of which the GLR equation could be a first ap- 
proximation. Again no scales are given and another of many hotly argued 
questions raised by the HERA data - for example in the context of the 
Caldwell plot (Fig. [7J) - is does HERA have the reach to see saturation 
effects and have they been seen? 

Dipole models have proved to be a very fruitful approach in exploring 
this question and modelling the behaviour of F2 or <r{^*p) through the 

4 An early version of this figure was shown by JK in his plenary talk on low x QCD at 
the 1993 Durham Workshop [H]. 
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transition region towards Q 2 = 0. The idea is that the virtual photon splits 
into a qq pair 5 (a colour dipole) with a transverse size r ~ 1/Q, the splitting 
occurring a 'long' time ~ l/{mx) before the dipole interacts with the proton. 
The dipole then scatters coherently from the proton in a time which is short 
in comparison to ~ l/(mx). The 7* — > qq process is described by QED and 
the strong interaction physics is then contained in the cross-section, a, for 
the dipole-proton interaction. 

Many authors have worked on these models in the context of DIS pro- 
cesses both inclusive and for diffraction. Here the discussion will focus on 
the model of Golec-Biernat and Wiisthoff (GBW) [Ej. The original GBW 
model provides a simple parameterization of a with explicit saturation, in 
the sense that a — > <jq (a constant) as r becomes large, but in a such a way 
that the approach to saturation is x dependent, controlled by the 'saturation 
radius', Rq. Explicitly the dipole cross-section is given by 

a(x,r 2 ) =a (l-exp(-r 2 )), f = — — (11) 

2R (x) 

where <ro is a constant, and 

Ro(x) 2 = ^(-)\ (12) 

Ro(x) 2 is understood to be inversely proportional to the gluon density (A ~ 
X g ), so that i?o is a measure of the transverse separation of the gluons in 
the target. Thus when the dipole separation is large compared to the gluon 
separation (small Q 2 and small x) the dipole cross-section saturates, a ~ cr , 
and <r(7*p) also tends to a constant, giving F2 cx 1/Q 2 from Eq.H When 
the dipole separation is small compared to the gluon separation (high Q 2 
and large x), a ~ <7o/(Q 2 -Ro) an d &{l*p) varies as 1/Q 2 , so that F2 exhibits 
Bjorken scaling. Because of the x dependence of Ro, a saturates for smaller 
dipoles sizes as x decreases. 

Looking at Fig. I1U1 which shows F2 as a function of Q 2 at fixed y = 
Q 2 /sx, one sees that the data do exhibit these features. There is a clear 
change in behaviour around Q 2 ~ 1 GeV 2 , which might then be taken as a 
rough estimate of the saturation scale for HERA data. However, one needs 
to be somewhat careful before jumping to the conclusion that the GBW 
dipole model shows that saturation effects have been seen at HERA. Given 
Eq.El F2 ~ 1/Q 2 at low Q 2 is required in any model since the photoproduc- 
tion cross-section is finite and slow, logarithmic scale breaking at high Q 2 is 
required by pQCD. Thus the real challenge is not in describing the general 



5 As for example shown by the quark box at the top of Fig. [3 
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Fig. 10. The GBW dipole model fit to F2 data through the transition region to 
photoproduction. The data are plotted versus Q 2 in bins of y. The dashed curves 
show the original model, the full curve the improved model with DGLAP evolution. 
From [Pd| . 



behaviour but in describing the exact shape of the data across the transi- 
tion, now that high precision data are available. It should also be noted 
that there are other dipole models that fit the same data as successfully 
without an explicit x dependent saturation scale. 

Although it gave a good description of F2 at low Q 2 , the original GBW 
model did not include QCD evolution in Q 2 and this limited its success at 
larger Q 2 . This deficiency was corrected in a second version |16| in which a 
was related to the gluon density. Results from both versions are shown in 
Fig. EH 
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Fig. 11. Geometrical scaling. Data on cr(^*p) with x < 0.01 plotted versus the 
scaling variable r = Q 2 R^{x). From [T7| . 

Another feature of saturation models which is illustrated in the GBW 
model is the dependence of the dipole cross-section on a scaled variable, 
here r/Ro(x). Stasto, Golec-Biernat and Kwiecinski ^7] have made the 
interesting observation that this leads to a new scaling property of a(j*p). 
At low x, cr(-f*p) should depend only on the dimensionless variable r = 
Q 2 Rq. What is more, such a scaling property seems to be satisfied by 
HERA data. One expects 

°~(l*p) ~ °~o ( T small) — > <r(7*p) ~ oo/r ( T large). 

Fig. ^2 shows HERA data with x < 0.01 as a function of r. Not only do 
the data show this scaling very clearly, but also the expected change in 
behaviour indicated above seems to occur around values of r ~ 1. This new 
scaling - geometrical scaling - holds for Q 2 values up to about 400 GeV 2 
but only for x < 0.01. Again, though it is not a proof, the remarkable 
demonstration of geometrical scaling by the low x HERA data does show 
that they have many of the attributes of a saturated system. 
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5. Summary 

The low x data from HERA have produced many surprising and interest- 
ing results that can be understood in the perturbative dynamics of a gluon 
rich system, although details such as exactly how to include BFKL effects 
are as yet undecided. Whether saturation effects have been demonstrated 
at HERA is controversial, but whatever the answer to this question one 
cannot deny the importance of the measurements in stimulating activity in 
this area. Particularly encouraging are the attempts to describe high den- 
sity gluon dynamics in terms of semi-classical models - the so-called 'colour 
glass condensate', which could be important for understanding data from 
heavy ion collisions at RHIC and the LHC. 

Jan Kwiecinski has been active in most of these fields and his work can 
be seen as an attempt to understand and set scales on the map shown in 
Fig. |§J It has been a pleasure to profit from his insights and untiring efforts 
to unravel the dynamic structure of the proton. 
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